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The pure Brans-Dicke (BD) gravity with or without the cosmological constant A has been taken 



2jJ! as a model theory for the dark matter. Indeed, there has been a consensus that unless one modifies 

■^ ' either the standard theory of gravity, namely, general relativity, or the standard model for particle 

pg ■ physics, or both, one can never achieve a satisfying understanding of the phenomena associated 

fT^ , with dark matter and dark energy. Along this line, our dark matter model in this work can 

lO ' be thought of as an attempt to modify the gravity side alone in the simplest fashion to achieve 

vn : 

^D ' the goal. Among others, it is demonstrated that our model theory can successfully predict the 

^^ . emergence of dark matter halo-like configuration in terms of a self-gravitating spacetime solution 

^^ I to the BD field equations and reproduce the flattened rotation curve in this dark halo-like object 

Qh' in terms of the non-trivial energy density of the BD scalar fleld, which was absent in the context 

5^ ■ of general relativity where Newton's constant is strictly a "constant" having no dynamics. Our 
-4— » ■ 

C^ , model theory, however, is not entirely without flaw, such as the prediction of relativistic jets in all 



types of galaxies which actually is not the case. 
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I. INTRODUCTION 

There has been a consensus among researchers that unless one modifies either the stan- 
dard theory of gravity, namely, general relativity, or standard particle physics theory (say, 
the Weinberg-Glashow-Salam's standard model), or both, one can never achieve a satisfying 
understanding of the phenomena associated with dark matter and dark energy. Along this 
line, in the present work, we would like to propose an attempt to modify the gravity side 
alone in the simplest fashion to achieve the goal. To be more specific, the pure Brans-Dicke 
(BD) gravity [l| with or without the cosmological constant A shall be taken as a model 
theory for the dark matter. Indeed, the BD theory is the most studied, and hence the 
3est-known, of all the alternative theories of classical gravity to Einstein's general relativity 



2|]. This theory can be thought of as a minimal 



extension of general relativity designed to 



properly accomodate both Mach's principle [2|, |3| and Dirac's large number hypothesis [2|, y]. 
Namely, the theory employs the viewpoint in which Newton's constant G is allowed to vary 
with space and time and can be written in terms of a scalar ("BD scalar") field as G = !/<!>. 
Besides, the BD scalar field (and the BD theory itself) is not of quantum origin. Rather, 
it is classical in nature and hence can be expected to serve as a very relevant candidate to 
play some role in the late-time evolution of the universe, such as the present epoch. 

As a scalar-tensor theory of gravity, the BD theory involves an adjustable, but unde- 
termined, "BD-parameter" u, and as is well-known, the larger the value of u, the more 
dominant the tensor (curvature) degree and the smaller the value of u, the larger the effect 
of the BD scalar. Also as long as we select a sufficiently large value of cj, the predictions of 
the theory agree perfectly with all the observations/experiments to date [j]. For this reason, 
the BD theory has remained a viable theory of classical gravity. However, no particularly 
overriding reason has ever emerged to take it seriously over general relativity. As shall be 
presented shortly in this work, here we emphasize that it is the expected existence of dark 
matter (and dark energy as well, see Ref.4) that may put the BD theory over general rela- 
tivity as a more relevant theory of classical gravity consistent with observations that have 
so far been unexplained within the context of general relativity. 

As mentioned above, in our model theory for dark matter, the effect of the cosmological 
constant A shall be generally considered. Here, A is essentially supposed to play the role 
of dark energy in the context of the BD theory as has been studied in detail in Ref.4. As 



is well-known, the mysterious flattened rotation curves observed for so long in the outer 
regions of galactic halos have been the primary cause that called for the existence of dark 
matter. Among others, therefore, we shall demonstrate in this work that our model theory 
can successfully predict the emergence of a dark matter halo-like configuration in terms 
of a self-gravitating static and nearly spherically-symmetric spacetime solution to the BD 
field equations and reproduce the flattened rotation curve in the outer region of this dark 
halo-like object in terms of the non-trivial energy density of the BD scalar field, which is 
absent in the context of general relativity where Newton's constant is strictly a "constant" 
having no dynamics. 



II. SCHWARZSCHILD-DE SITTER-TYPE SOLUTION IN THE BD THEORY OF 
GRAVITY 

As stated earlier in the introduction, we would like to demonstrate in the present 
work that the BD gravity with or without the cosmological constant can reproduce some 
representative features of dark matter, such as the formation of a dark matter halo 
inside of which the flattened rotation curves are observed. Since the galactic dark matter 
halos are roughly static and spherically-symmetric, we should, among others, look for 
such dark matter halo-like solution to the BD field equations. Thus, in this section, 
we shall construct the Schwarzschild-de Sitter-type solution in the BD theory and claim 
later on that it can represent the dark matter halo in the context of our model of dark matter. 

As is well-known, even in Einstein's general relativity. Finding exact solutions to the 
highly non-linear Einstein field equations is a formidable task. For this reason, algorithms 
generating exact, new solutions from the known solutions of simpler situations have been 
actively looked for, and actually quite a few have been found. In the BD theory of gravity, 
the field equations are even more complex; thus, it is natural to seek similar algorithms 
generating exact solutions from the already known simpler solutions either of the BD 
theory or of the conventional Einstein gravity. In particular, Tiwari and Nayak [5| proposed 
an algorithm that allows stationary axisymmetric solutions in the vacuum BD theory to 
be generated from the known Kerr solution in the vacuum Einstein theory. Thus, in the 



present work, we shall take the algorithm suggested by Tiwari and Nayak [5] or by Singh 
and Rai |6| to construct the Schwarzchild-de Sitter-type solution in the BD theory in the 
presence of the cosmological constant from the well-known Schwarzschild-de Sitter solution 
in Einstein gravity with the cosmological constant. 



Consider the BD theory in the presence of the (positive) cosmological constant A de- 
scribed by the action 



S 



d^Xy/g 



1 /^„ V„$V°<I> 

^R-LU 



A 



(1) 



.167r V $ 

where $ is the BD scalar field representing roughly the inverse of Newton's constant and u 
is a generic parameter of the theory. Extremizing this action with respect to the metric Qf^iy 
and the BD scalar field $ yields the classical field equations given, respectively, by 
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The Einstein gravity is the u; ^ oo limit of this BD theory. Now the algorithm of Tiwari 
and Nayak [5| or Singh and Rai [6] goes as follows; In general, let the metric for a stationary 



axisymmetric solution to the Einstein field equations take the form 
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while letting the metric for a stationary axisymmetric solution to the BD field equations be 
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where U, W, k, and h are functions of x^ and x^ only. The significance of the choice of the 
metric in this form has been thoroughly discussed by Matzner and Misner [7| and Misra and 
Pandey |8|. Tiwari and Nayak or Singh and Rai first wrote down the Einstein and the BD 
field equations for the choice of metrics in Eqs. (3) and (4), respectively. Comparing the 
two sets of field equations closely, then, they realized that stationary axisymmetric solutions 
of the BD field equations are obtainable from those of the Einstein field equations provided 
certain relations between metric functions hold. That is, if {We, kE, Ue, He) form a 
stationary axisymmetric solution to the Einstein field equations for the metric in Eq. (3), 



then a corresponding stationary axisymmetric solution to the BD field equations for the 
metric in Eq. (4) is given by {Wbd, ksn, Ubd, hsD, $), where 
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$ = [/i^]4/(2-+3). 



Now what remains is to apply this method to obtain the Schwarzschild-de Sitter-type 
solution in the BD theory from the known Schwarzschild-de Sitter solution in Einstein 
gravity. To do so, one needs some preparation, which involves casting the Schwarzschild-de 
Sitter solution given in the usual Schwarzschild coordinates (t, r, 6, (p) in the metric form 
given in Eq. (3) by performing the coordinate transformation (of r alone) suggested by 
Misra and Pandey 81]. Namely, we start with 
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where A = r^(l — 87rAr^/3) — 2Mr with M and A denoting the ADM mass and the (posi- 
tive) cosmological constant, respectively. Consider now such a transformation of the radial 
coordinate as the one introduced by Misra and Pandey [8|; 

dr'^ 
r = r(R) such that — — = dR^. 
^ ' A 

Then, the Schwarzschild-de Sitter solution can now be cast in the form in Eq. (3), i.e., 
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with L = r{R) and, hence, A = L^(l — 87rAL^/3) — 2ML. Now, we can read off the metric 
components as 
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Then, using the rule in Eq. (5) in the algorithm by Tiwari and Nayak or by Singh and Rai, 
one finds that the metric components of the Schwarzschild-de Sitter-type solution in the BD 
theory come out to be 
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Next, looking at the forms of the metric for a stationary axisymmetric solution to the 
Einstein or the BD field equations given in e qs.( 3) and (4), one can realize that the factors 
in the metric functions, 
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should have length dimension 2. Firstly, for the solution to the Einstein field equations. 
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which all take the right length dimension, 2, as expected. Secondly, for the solution to the 
BD field equations, however, we have 
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Thus, in order for this metric solution to the BD field equations to have the right 
length dimension, it turns out that the factor (Asin2 ^)'^2/(2u;+3) gi^ould be dimension- 
less. Therefore, in order to guarantee this, we introduce a normalization factor tq to 
render the factor A (particularly appearing in (Asin^ ^)^2/(2i^+3)^ dimensionless, i.e., 
A = [L'^{1 — 87rAL2/3) — 2ML]/rQ. The proper value of this normalization factor Tq shall 
be characterized later on when we do the numerics to quantify the rotation velocity in our 
model. 

Then, by transforming back to the standard Schwarzschild coordinates using Eq. (7), we 
finally arrive at the Schwarzschild- de Sitter-type solution in the Brans-Dicke theory given 
by 
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where we have restored both the present value of Newton's constant Gq and the speed 
of light c in order to come from the geometrized unit (Go = c = 1) back to CGS 



6 



units. That is, we have replaced M with G^Mj(? = M along with AGqA/c^ = A;hence, 
A = [r^(l — 87rAr^/3) — 2Mr]/r'^. A remarkable feature of this Schwarzschild-de Sitter-type 
solution is the fact that, unlike the Schwarzschild-de Sitter solution in general relativity, 
the spacetime it describes is static (i.e., non-rotating), but not spherically-symmetric. Note 
also that as cj -^ oo, this Schwarzschild-de Sitter- type solution goes over to the standard 
Schwarzschild-de Sitter solution in Einstein gravity as it should since the uj ^ oo limit of 
the BD theory is the Einstein gravity. Note also that the Brans-Dicke-Schwarzschild (EDS) 
spacetime solution to these vacuum BD field equations has been found some time ago |9| in 
rather a theoretical attempt to construct non-trivial black hole spacetime solutions in the 
BD theory. There the EDS spacetime solution has been obtained by setting a = e = in 
the Brans-Dicke-Kerr-Newman (BDKN) solution in Eq. (11) of Ref.9, and it coincides with 
the Schwarzschild-type solution that results when we set A = in the solution constructed 
above, as it should. 

Now coming back to our main objective, we are interested in the role played by the 
BD scalar field as dark matter, particularly in forming galactic dark matter halos inside 
of which the well-known rotation curves have been observed. We suggest (and later on 
demonstrate in terms of the reproduction of a flat rotation curve) that, indeed, the BD 
scalar field can successfully play the role of dark matter. In order for this to happen, 
we first claim that the BD scalar can cluster into a halo-like configuration as it can be 
represented by this Schwarzschild-de Sitter-type solution constructed above in Eq. (14). 
At this point, therefore, it seems relevant to address the nature of potential singularities 
of this Schwarzschild-de Sitter-type metric solution. Just as the Schwarzschild-de Sitter 
solution in general relativity, it appears to possess two coordinate singularities which would 
arise at A = [r^(l — SvrAr^/S) — 2Mr]/rQ = ; (i) the inner Schwarzschild gravitational 
radius, and (ii) the outer de Sitter radius. Generally, it is well-known that the Schwarzschild 
metric solution in static coordinates possesses an event horizon, that the solution essentially 
describes the region outside this event horizon, and that its inside is not everywhere well- 
defined. Meanwhile, the de Sitter metric solution in static coordinates has a cosmological 
horizon,the solution mainly describes the region inside this cosmological horizon, and its 
outside is not everywhere well-defined. As such, the Schwarzschild-de Sitter metric solution 
in static coordinates is able to represent the regionbetween the inner event horizon and the 



outer cosmological horizon, but not elsewhere. In the same spirit, our Schwarzschild-de 
Sitter-type metric solution in the BD theory given in static coordinates has been adopted 
here only to represent the galactic dark halo region, but not below (i.e., the interior of 
a given galaxy) nor beyond (i.e., the scale of galaxy clusters or even of entire universe). 
Therefore, all we have to do is to demonstrate that if we quantify the locations of these two 
singularities by putting real numbers, they are well below and beyond the typical domain 
of galactic halos. 

First, we start with (i) the inner Schwarzschild gravitational radius. For small-r, A = 
is approximated by r — 2M ~ whose solution is r ~ 2GqM/(:? ~ 0.01 pc. Here we used, 
for the typical (total) mass of a galaxy, M ~ lO^^M©. Obviously, this occurs well inside 
a galaxy as the typical size of a galaxy ranges from few kpc (for dwarf galaxies) to few 
hundred kpc (for ordinary galaxies). Indeed, the physical meaning of this singularity is that 
if the entire galaxy is squeezed into this gravitational radius, it becomes a black hole with 
its event horizon placed at 0.01 pc. Namely, for ordinary galaxies, this inner Schwarzschild 
gravitational radius is a failed Schwarzschild event horizon. Next, we turn to (ii) the outer 
de Sitter radius. For large-r, A = is approximated by r — SvrAr'^/S ~ 0, whose solution 
is r ~ (SttCqA/Sc^) ^ 4 Gpc. Here we used for the cosmological constant the observed 
value, A ~ 10"^ ergjcw? . Apparently, this occurs well beyond a single galaxy; Indeed, this 
is the scale of entire universe (given the age of the universe, that is r ~ 13.7 Gyr, its rough 
size would be cr ~ 4 Gpc). In other words, this outer de Sitter radius is totally irrelevant. 

To summarize, the potential singularities of the Schwarzschild-de Sitter-type solution are 
irrelevant to keeping us from employing the solution to represent the galactic dark halo region 
and, hence, are harmless. Next, the seeming angular singularity at ^ = 0, vr are irrelevant 
as well because the domain of principal interest is the neighborhood of galactic equatorial 
plane, 6 = 7r/2, inside the halo where most gases and stars orbit around with flat rotation 
curves. Normally, the symmetry axis is the last thing to be expected to be singular. Indeed, 
close inspection reveals that the metric function factor (Asin^6')^^/''^'^"^^'* is generated by 
the solution-generating algorithm of Tiwari and Nayak pj or of Singh and Rai 61] which 
takes the (singularity-free) solutions of the Einstein field equations to those of the BD field 
equations. That is, since the appearance of the factor (sin^^)^^/*^^'^+^\ which is responsible 
for the singular nature of the symmetry axis, results simply from the solution-generating 
algorithm, one naturally might expect that the symmetry axis ^ = 0, vr cannot possibly be a 
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genuine curvature singularity as it would not really represent, say, an infinite concentration 
of matter along there. However, the symmetry axis ^ = 0, vr, indeed, appears to be singular 
as the invariant curvature polynomial, such as the curvature scalar R, which can be readily 
calculated by contracting the metric field equation and using the field equation for the BD 
scalar field in Eq. (2), is given by 

R = 7T^^4r(^ «i^' ^)^2/(2.+3) [2^.2 + 4 1(2^ + 3)AnAGy + cot^ O} A , 

where A' = 2r 1 — (M/r) — (IGvrAr^/S) /r^, blows up there for generic u values. Namely, 
the symmetry axis appears to be a real singularity rather than a coordinate singularity 
(that can be removed by a change of coordinates). This certainly is a very unexpected and 
hence puzzling feature of the Schwarzschild-de Sitter-type solution in BD gravity theory. 
Therefore, one might wonder what would happen to matter habiting near the galactic poles 
when it comes close to the pole. We now hope to clarify the true nature of this peculiar 
singularity along the symmetry axis in some detail. 

Among others, the first thing that comes to our mind regarding the nature of this sin- 
gularity along the symmetry axis is the fact that, unlike the familiar curvature singularity 
at the center of a black hole, which is point-like, this singularity is an infinitely-extended 
line singularity. This means that perhaps the fate of matter coming close to this singularity 
would be somewhat different from what one would normally expect for matter approaching a 
point-like black hole singularity. To get right to the point, it turns out that in the immediate 
vicinity of the symmetry axis, the specific energy (i.e., the energy per unit mass) of a test 
particle becomes extremely high. As a result, the particle moves along the symmetry axis at 
nearly the ultrarelativistic speed (i.e., at nearly the speed of light). Of course, this is an un- 
usual feature, which has no analogue in the Einstein gravity context, and a close inspection 
reveals that it can eventually be attributed to the metric function factor (A sin^ ^)±2/(2<^+3)^ 
which is responsible for the singular nature of the symmetry axis. The explicit demonstra- 
tion of this peculiar behavior of test particles near the symmetry axis shall be presented 
later on in section V in terms of a rigorous analysis of the geodesic motion there. 

In conclusion, this study of the true nature of the singularity along the symmetry axis 
leads us to suspect that perhaps the bizzare singularity at ^ = 0, tt of the Schwarzschild-de 
Sitter-type solution in BD gravity theory can account for the relativistic bipolar outflows 
(twin jets) extending from the central region of "active galactic nuclei (AGNs)" Namely, 



the curious singularity along the symmetry axis seems harmless, after all. Instead, it turns 
out to be a pleasant surprise as it can explain a long-known puzzle in observed features 
of galaxies. A cautious comment might be needed here, though. That is, the relativistic 
bipolar outflows have been observed only for some particular types of galaxies, such as 
AGNs and micro- quasars. (Radio galaxies, Seyfert galaxies and quasars fall into the AGN 
category.) Namely, the jets do not seem to be a general feature of all types of galaxies. 

Lastly, the failure of asymptotic flatness of this Schwarzschild-de Sitter-type solution is 
not of serious concern here as we essentially aim at flnding a spacetime solution that can 
represent the dark halo-like configuration, which is known to cluster only on a snh-Mpc 
scale, i.e., a local scale, and eventually matches the cosmological geometry, such as the 
Friedmann-Robertson- Walker metric for homogeneous and isotropic expanding universe, at 
a larger cosmological scale in order for the BD theory in the presence of the cosmological 
constant to provide a successful model for dark matter and dark energy, as suggested, for 
example, in Ref.4. 

Our natural next mission is then to ask whether these configurations really can reproduce 
the properties of dark matter halos, namely, if our BD scalar model for dark matter can 
reproduce the fiattening of the rotation velocity curves inside these halo configurations 
consistent with the observations. Therefore, in the following, we shall address this issue; 
first, for the case when only the dark matter content represented by the BD scalar is present, 
i.e., when the dark matter halo is represented by the Schwarzschild-type solution (without 
the cosmological constant (A) term), andthen for the other case when both the dark matter 
content and the dark energy content, represented by the cosmological constant, are present, 
i.e., when the dark halo is represented by the Schwarzschild-de Sitter-type solution (with 
the A term). 

III. BD DARK MATTER HALOS NEGLECTING THE DARK ENERGY (A) 
CONTRIBUTION: THE CASE OF THE JORDAN FRAME 

As claimed earlier, it appears that the BD scalar can indeed cluster into a halo-like 
configuration as it can be represented by the Schwarzschild-type solution in the BD theory. 
Thus, we now attempt to obtain the rotation curves in our BD scalar halo and eventually to 
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demonstrate that they are actually flattened far out to the distant region of halo. First of all, 
since we need concrete "numbers" we now restore both Newton's constant Go and the speed 
of light c in order to come from the geometrized unit (Gq = c = 1) back to the CGS unit. 
Then, the energy-momentum tensor of the BD scalar field given earlier in Eq. (2) should 
now be multiplied by the factor (c^/Gq), and in the Schwarzschild-type solution in eEq. 
(14) above, we should replace M by GqM/(? = M, as we mentioned earlier. Apparently 
then, the energy-momentum tensor of the BD scalar field with restored Go and c has the 
dimension of the energy-momentum density in CGS units {erg /civ?). 

We now turn to the computation of energy density profile and (anisotropic) pressure 
components of the BD scalar field playing the role of the dark matter by treating the BD 
scalar field as a (dark matter) fluid. The BD scalar field fiuid, however, would fail to be a 
"perfect" fiuid as can readily be envisaged from the fact that the associated BDS solution 
configuration is not spherically-symmetric. Namely, its pressure cannot be "isotropic", i.e., 
P,. ^ Pe ^ P(f,. Such a fiuid may be called an imperfect fiuid due to the anisotropic pressure 
components, and as such, its stress tensor can be written as 

/-cV \ 
Pr T; 
T^^ Pe 

V p^y 

which is to be contrasted to its counterpart for the usual perfect fiuid with isotropic pressure 
given by the well-known form Tj^ = P6^ + {c^p + P)U^Uu = diag{—c^p, P, P, P), where 
f/" = dX°'/dT (with r being the proper time) denotes the 4- velocity of the fiuid element 
normalized such that If^Ua = —(? . Note that in addition to the diagonal entries representing 
the (anisotropic) pressure components Tl = Pi. (with no sum over i), there are off-diagonal 
entries Tg, Tj! representing a shear stress that also results from the failure of spherical 
symmetry. It is also interesting to note that a stress tensor of this sort (given in eEq. (15)) 



n^BD tj. 



(15) 



arises in the case of rotating boson star made up of a complex scalar field ll|. Thus by 
substituting the BDS solution given in Eq. (14) into the BD energy-momentumm tensor in 
eq.(2) and then setting it equal to Eq.(15), we can eventually read off the energy density 
and the pressure components of the BD scalar field imperfect fiuid; 
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where A = [r{r — 2M)]/rQ and A = TqA. Note that the off-diagonal components Tg and 
T^ are odd functions of 9 while the diagonal components (p, P^, Pe, P4) are even functions 
of the polar angle under 9 ^ {-k — 9). As a result, the off-diagonal components vanish (i.e., 
no shear stress survives) if we average over this polar angle to get a net stress. Thus, the 
equation of state of this BD scalar k-essence fluid forming a galactic halo is given by 
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c^p 
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2{uj + 1) {(r - M)2 + A cot2 9] - (2u + 3)M(r - M) ' 
where P = Pr- Namely, P = w{r,9)c^p with w{r,9) ~ 0(1), meaning that this BD scalar 
fluid is essentially a barotropic fluid but with a "position-dependent" coefficient w{r,9). 
Note that although the BD scalar field is a candidate for dark matter, it is not quite a 
dust. In principle, the speed of sound in this BD scalar field fluid can also be evaluated via 
c2 = dP/dp, but we shall not discuss that in any more detail in this work. 

We are now ready to compute the behavior of the rotation curves in the outer region (i.e., 
at large, but finiter, say, r >> GqM/c^ ) of our BD scalar field halo. To be more precise, 
for a galaxy of typical (total) mass M ~ IO^^Mq, the outer region of its dark matter halo, 
say, r ~ 10 (kpc) ~ 102'^ (cm) is much greater than GqM/c^ ~ 10^^ {cm) by a factor of 10''' 
or so. Thus, to this end, we first approximate the expressions for its energy density and the 
(radial) pressure given in Eq. (16) for large-r; 
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Note that in the above approximations and in the discussions below, it was and it shall be 
assumed that the metric function A = [r{r — 2M)]/rQ ~ {t/tqY for larger. It is interesting 
to note that as a candidate for dark matter, the energy density p of the BD scalar field is 
almost certainly positive everywhere (i.e., for both small and larger). In the meantime, its 
(radial) pressure P particularly on a larger scale (i.e., for larger) turns out to be negative 
although its sign appears unclear on a small scale (i.e., for smallr). 

Finally, we are ready to determine the rotation curve inside our BD scalar field halo. 
First, we start by recalling the origin of the long-standing puzzle associated with the galaxy 
rotation curves. In the most naive sense, the apparent rotation velocity of an object at a 
radius r from the galactic center would be given in the Newtonian limit by Kepler's third 
law, f2 = GoM{r)/r. Indeed, inside a given galaxy where the luminous mass is nearly 
uniformly distributed, the rotation velocity has been observed to grow roughly linearly with 
the distance r, consistent with this Kepler's law. In the outer region of a galactic halo, 
however, the rotation velocity is expected to behave like v ~ 1/v^ ^^ the luminous mass is 
confined within the extent of the given galaxy. The observations, however, exhibit flattened 
rotation curves, and this has been the age-old dilemma that called for the existence of dark 
matter in the outer regions of galactic halos. 

In the present work, therefore, we shall employ the expression for the rotation velocity 
still given by Kepler's law, but instead attempt to explain the flattened rotation curve in 
terms of the non-trivial energy density of the BD scalar field, which appears to play the 
role of dark matter. Namely, the non-trivial mass density of the BD scalar field turns out 
to contribute to the mass function M{r) in such a way as to render the rotaion curve flat, 
as we shall see in a moment. 

At this point, regarding the construction of the rotation velocity, we have some 
comments. For the present case, it may seem that one cannot just employ Kepler's law 
as the stress tensor of the BD scalar fleld (playing the role of dark matter) and that the 
background Schwarzschild-de Sitter-type spacetime in the BD theory is not spherically 
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symmetric. Nevertheless, gases or stars orbiting around the host galaxy center usually lie 
on the galactic equatorial plane. This implies that the orbital motions on the equatorial 
plane 9 = it/2 are of particular interest; hence, there one may still wish to employ Kepler's 
law. Along this line, therefore, one might wish to employ the relativistic counterpart 
to this rotation velocity equation, instead. Later on in the appendix, the construction 
of the rotation velocity shall be promoted to a fully relativistic version in terms of the 
rigorous derivation of timelike geodesies in the Schwarzschild-de Sitter-type spacetime. 
The quantitative results that we are about to present below, however, essentially remain 
unchanged even for the fully relativistic treatment. 

Now for our case, using the BD scalar energy density profile 
given earlier, we have M{r) = Jq"" dcj) J^~' dO /J" dr^g„geeg<f>4>p{r, 0) = 
{2^/Gq) \{uj + l)/{2uj + l){2uj + 3)] /(c.;)r(r/ro) -2/(2^+3) and hence, 

2 , , _ GQM{r) _ ^2 2(a; + 1) ^^ .^ ( r ^-2/(2-+3) 



^ ' r (2cj + l)(2cu + 3)^^ '\ro) ' ^ ^ 

where f{uj) = /^"^ rf^sin-[^+2/(2-+3)l = 2 J^'^ dx[l-x^]-^^'^+^^/^^^+^^ withe, S « 1. (Note 
here that the integration over the polar angle 6 starts not from but from e << 1 as the 
symmetry axis 6* = of the BDS solution in Eq. (14) possesses the danger of an internal 
infinity nature, namely, the symmetry axis is an infinite proper distance away, as discussed 
carefully in Ref.lO.) It has been known for some time that in order for the BD theory to 
remain a viable theory of classical gravity passing all the observational/experimental tests to 
date, the BD cu-parameter has to have a large value, say, \uj\ > 500 [2(]. In our previous study 
9j, however, we realized that the static solution to the vacuum BD field equations given in 
Eq. (14) above, but without A, can turn into a black hole spacetime for —5/2 < u < —3/2. 
Thus, now for |ti;| > 500, the same static solution Eq. (14) we are considering represents just 
a halo-like configuration with a regular geometry everywhere (i.e., having no horizon), which 
is static but not exactly spherically symmetric (note that the galactic halos are also believed 
to be nearly spherically symmetric, but not exactly). Thus, if we substitute a large-u; value, 
say, uj ~ 10^ into Eq. (19) above, evidently M{r) ~ r; hence, we get 

/ ^ s -(1/106) 

v{r) ^ 100(A;m/s) x — (20) 

Vro/ 

because for u ~ 10^, f{uj) — 0(1). Namely, for this large-cj value, the rotation curve gets 
flattened out as r~(^° ^ ~ constant, and its magnitude becomes several hundred km/ s. 
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Here, it is particularly remarkable that the observationally and experimentally allowed 
large value of tu ~ 10^ renders the rotation curve flat and fixes its magnitude to several 
hundred km/ s at the same time! This is, indeed, an attractive feature of the BD theory 
that distinguishes it from other theories of dark matter proposed thus far. Of course, this is 
in impressive agreement with the data for rotation curves observed in spiral/elliptic galaxies 
with M/ L ~ (10 — 20)Mq/Lq and in low-surface-brightness (LSB)/dwarf galaxies with 
M/L ~ (200 — 6OO)M0/L0 (where M/L denotes the so-called "mass-to-light" ratio given in 
the units of the solar mass-to-luminosity ratio and exhibiting a large excess of dark matter 
over the luminous matter) [12]. For instance, the rotation curve of the dwarf spiral galaxy 
M33 is shown in Fig.l. Rotation curves are observed usually via measurements of the 
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FIG. 1: Flattened rotation curve in the dark halo of dwarf spiral galaxy M33. 

Doppler shift of the 21 — cm emission line from neutral hydrogen (HI) for distant galaxies 
and of the light emitted by stars of nearby galaxies |l3l . Il4| . It is also interesting to note that 
this behavior of the rotation curve in our BD theory dark matter halo model is independent 
of the mass of the host galaxy. Namely, this behavior of the rotation curve comes exclusively 
from the nature of the dark matter, i.e., the BD scalar field. Indeed, we have restricted 
our interest in the present work only to the contribution coming from the dark matter 
component, i.e., the BD scalar field and not to the combined contributions from both the 
luminous and the dark matter. Therefore, the resulting galactic rotation curve (having a 
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contribution coming only from the BD scalar, which is the dark matter candidate) comes 
out as being independent of the luminous mass in each galaxy. The well-known Tully-Fisher 



relation [15|] between the total luminosity of a given galaxy and the outermost rotation 
velocity can be accounted for if both the luminous and the dark matter contributions are 
considered. This issue has not been addressed in the present work, but hopefully it should 
be answered by any dark matter model in order for the proposed model to be truly successful. 

We also point out that even if we employ more exact expression for the rotation velocity 

curve involving the Doppler shift of light emitted by the orbiting objects (assuming that the 

^-essence halo is almost sphericallysymmetric) , namel,y v'^{r) = GQM{r)/r + Anr'^GoP/c^ 

16| (with P being the radial pressure given in Eq. (18) above), the conclusions above remain 

the same. Thus, in what follows, we shall demonstrate this in some detail. Using Eqs. (18) 

and (19), we have 

2 GoMjr) , AttGq ^ 
V (r) = 1 —r P (21) 

_ „. 2(J+ 1) \ ( r N -V(^-+3) ^^ 2[2(a; - 1) cot^ 6 + 1] .. ■4/(2.+3) , ( r y'^^'-^'^ 

Again, due to the failure of spherical symmetry in the BDS solution, Eq.(14), which is 
identified with a dark matter halo configuration, this rotation velocity comes out to be polar 
angle (9) dependent. However since the gases or the stars orbiting around the galaxy center 
usually lie on the galactic (equatorial) plane, the relevant situation to consider is the case 
9 = -n 12 in which 



f (r) 



2(cc; + 1) ( r\ -2/(2^^+3) 2 ( r\ -4/(2^^+3) 



r/(^) - 



(2cj + l)(2tu + 3) ' ' Vro/ (2cJ + 3)2Vro 



1/2 



r \-(i/io«) 
100(fcm/s) X ( — , (22) 



where again we have chosen the value lo ~ 10® in the last line. 

Next, the equation of state in eEq. (17) of this BD scalar field fluid becomes, in the outer 
region of the galactic dark matter halo (i.e., at larger), 

2(cj-l)cos2^ + sin2^" 



w 



(23) 



2(cj + l) 

which is obviously negative due to the negative pressure (and still positive energy density) 
in this outer region. Moreover, for a large-a; value, i.e., uj ~ 10®, for which the rotation 

16 



curve gets flattened out, that we just have reahzed, this equation of state at larger further 
approaches lu ~ — cos^ 9 c^ —0(1). (Incidentally, it is interesting to note that in the vicinity 
of the equatorial plane 6 = 7r/2, w = 0; namely, the BD scalar behaves nearly as dust.) This 
observation is particularly interesting as it appears to indicate that the BD scalar fleld we are 
considering possesses dark -energy-like negative pressure on larger scales. This observation is 
indeed consistent with our previous study (j|, which showed that on a cosmological scale, the 
BD scalar field, which has been identified with a "k-essence" there, did exhibit the nature 
of dark energy possessing a negative pressure. 

IV. BD DARK MATTER HALOS INCLUDING THE DARK ENERGY (A) CON- 
TRIBUTION 



By substituting the Schwarzschild-de Sitter-type solution given in Eq. (14),, but this 
time into the total energy-momentum tensor including the contribution coming from the 
cosmological constant term 
1 



T, 



flU 



T. 



BD 



Gn^ 



A^, 



flU 



(24) 



SttGo 



|^(V^$V,$ - ^(?^,V„$V"$) + ^(V^V.'I' - ^^.V^V'^'I') 



1 



Go^ 



A^, 



fit^i 



and then again setting it equal to Eq.(15), we can eventually read off the energy density and 
the pressure components of the BD scalar field imperfect fluid as 



P 



1 /a • 2/,\-2/(2-+3) 

A sm 6 



{r- M) Ar^ 



SttGo (2cj + 3)2 r2A 
2uj {{r-M) -ArA + (2(r - M) - (2cj + 3)M + {2u - 1)— Ar^ 



+ 2(cj + l)Acot^^ 
c* 4 



+ 



(^ 



sin^^ 



X 



SttGo (2cj + 3)2 r2A 
2(u; + l)((r-M)- 



(2tu - 1) A . . 2 
(2a; + 3) c^ 

1 U ■ 2 /,\ -2/(2^+3) 

A sm^ 9 



-2/(2a;+3) 



r - M) Ar^ 



— Ar4 + |2(r - M) - (2cu + 3)M + {2u - 1)— Ar=^ 



- (2cj + 3)A(l -167rAr^) + 2(cj-l)Acot2 



(2w - 1) 
(2w + 3) 



A (a sin' Q) 



-2/(2w+3) 



SttGo (2w + 3)2 r2A 
{r -M) Ar^ 



LfAsin'e 



-2/(2w+3) 



(2cj + 3) 2{uj - 1) |(r - M) ^Ar^ 
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{2(cj+l)cos^^-(2cj + 3 



A 



sin^^ 



{2lu - 1) 
(2cj + 3)' 



A Asin^^ 



-2/(2w+3) 



(25) 



SvrGo (2cj + 3)2 r2A 



-2/(2a;+3) 



X 



TJ 



(r-M) 
A cot^ 9 

AT' 



3 

(2c. - r 



A sin^ 9 
2{uj + 1) j (r - M) 

2/{2a;+3) 






2cu + 3 — 
r 



(2cj + 3) 

:4 4 



AfAsin^^ 
1 



SttGo (2cu + 3)2 r2 
4cj <! (r - M) 



cot 6* A sin^ 9 



-2/(2w+3) 



i^Ar4-(2c. + 3)^ 
3 J r 

where now A = [r^{l — 8'KAr^/3)—2Mr]/rQ and A = TqA. As in the previous case where the 
contribution from the dark energy component, namely, A has been neglected, the off-diagonal 
components Tg , T^! are odd functions of 9 whereas the diagonal components {p,Pr,P0,P(f,) 
are even functions of the polar angle under 9 -^ (vr — ^). Consequently, the off-diagonal 
components vanish (i.e., no shear stress survives) if we average over this polar angle to get 
a net stress. Again, first the equation of state of this BD scalar k-essence fluid forming a 
galactic halo is given by 



w 



c^p 



(26) 



Equations (25) and (26) indicate that P = w{r, 9)c^p with w{r, 9) ~ 0(1) meaning that this 
BD scalar fluid is still a barotropic fluid but with "position-dependent" coefficient w{r,9). 

As before, we now turn to the exploration of the behavior of rotation curves in the outer 
region (i.e., at large but finite-r, say, r ~ 10 (kpc) ~ 102^ (cm) » GqM/c^ ~ 10"*^^ {cm)) 
of our BD scalar field halo. And to this end, we first approximate the expressions for 
the energy density and the (radial) pressure of the BD scalar fluid given in eq.(25) using 
essentially A = [r2(l — 87rAr2/3) — 2Mr]/rl ~ [r2(l — 87rAr2/3)]/rQ for large-r. They are 



[UJ 



1) 



+ 



SttGo {2uj + 3)2 r2 sin^ 9 

{2uj - 1) 4 (2cj + 1) 



A sin^ 9 



-2/(2w+3) 



P 



{2uj + 3) 3 {2lu + 3)2_ 

C^ 1 1 r 

' 27iGo {2u + 3)^^ ^ ^'^ 
(2cj - 1) 4 1 



?(- 



sm 



-2/(2w+3) 



(27) 



l)cot2^+ 1 



A sin^ 9 



-2/(2ui+3) 



{2uj + 3) 3 (2u; + 3)^ 



A (a sin^ 9 



-2/(2a;+3) 



18 



It is interesting to note that the energy density p of the BD scalar field even in the presence 
of A, the dark energy component, is almost certainly positive everywhere (i.e., for both small 
and large-r). In the mean time, its (radial) pressure P, in the presence of A, particularly at 
larger scale (i.e., for large-r) turns out to be negative. 

We are now ready to determine the rotation curve inside our BD scalar field halo. First, 
using the BD scalar field energy density profile given earlier, we are supposed to com- 
pute the mass function M{r) = /q ^ dcj) J^~^ dO /J" dr^grrgeeQcfxfiPi'^^ &) in which we have 
y/9rr9ee94>4> = (Asin^6'j A^-*^/^ (r^/rg) sin ^. Although we shall consistently employ 

the approximation A ~ [r^(l — 87rAr^/3)]/rQ as we are interested in the behavior of rotation 
curves in the outer region (i.e., at large-r) of the BD scalar halo, the actual computation of 
the mass function M{r) is unfortunately not available since the analytic integration over r 
cannot be done in a closed form. Indeed one of our main objectives in this section, where 
the contribution from the dark energy component (i.e.. A) has been included, involves the 
theoretical derivation of, say, the upper bound on the value of the cosmological constant 
A that would result in, despite its presence, the nearly fiattened rotation curves as we ob- 
serve them. Obviously, the (positive) cosmological constant has positive contribution to the 
mass-energy density p and as a result, in the presence of A, the mass function M{r) grows 
non-linearly with the distance (i.e., the radius of the dark halo) ruining the asymptotic 
flatness of the rotation curves. As such, in order for our model, i.e., the BD theory with 
(positive) cosmological constant to successfully reproduce the flattened rotation curves, it is 
evident that the absolute magnitude of the cosmological constant should be small enough. 
For this reason, we shall, for the sake of explicit, analytic computation of the mass function 
M{r), assume that the dimensionless quantity, Ar^ be very small, i.e., (G'oA/c"^)r^ << 1 and 
see what we would end up with. Then we have 



M{r) 



2lT 







C 



i^+l] 



r 



dO / dr 

Jo [8716-0(20; + 3)2 sin ^ 



(2cu-l) 4 (2cu+l) 



(2cu + 3) 3 {2uo + 3)2 



^ 3 . / 

— r sm U 
c^ 



X A-i/2(Asin2^ 
c2 2iu + 1 



-l/(2w+3) 



Go {2uj + 3)2 



/(^) 



(2^ 



X 



271 



{2uj + 3) {2uo + 5) 47rGo 
(2cu+l)^^ (6cj + 7) 3c4 
{2uj - 1) 4(2cj + 1) I 
{2uj + 3) ^ 3(2cu + 3)2f 6 



Ar'^ + 



;iOcj + 13) V 3c' 



AtiGv. 



'-A 



-2/(2w+3) 



-9{^) 
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X 



(2a; + 3) 



r-" + 



(2cu + 5) 47rGc 



Ar- 



MttGo 



A r' 



J, \^ -2/(2a;+3) 



{Qlo + 7) {IOlo + 13) 3c4 (14cj + 19) V 3c^ 

where now /(cj) = J^-' desm-^^+^/^^'^+^^^ 9 = 2 J^~Ux[l - a;2]-(2^+4)/(2^+3) ^ 
0(1) and ^(cj) = /;~^rf^sin[i-2/(2..+3)]^ _ 2/o'-'^rfa;[l - x^]'^/'^^''+^^ ~ 
0(1) with e, 6 « 1. And here we used A^^/^^^^/^^o^+s) ^ 

(l + 47rGoArV3c4)(H-87rGoArV3(2cu + 3)c^)r-i(r/ro)-2/(^'^+^) since we assumed 
(OoA/c'^)r^ << 1 as explained above. Lastly, the rotation velocity is given by 
GoM(r) 2 2(cu + l) / ^ X -2/(2..+3) 



V (r) 



(2o; + l)(2tu + 3)^^'^- Vro 



+ 



+ c^ 



4(2cu+l) 



(2o; - 1) 

(6cj + 7) ~ 3{2uj + 3){6lu + 7) 

4(cj + 1) MttGo 

(2cu + 3)2(10cj+13) V 3c4 



27rOr 



Ag{uj) 



2(cj+l)(2cj + 5) 47rGo 
(2cu + 3)2(6cu + 7) 3c2 



(29) 



A/(c.) r' 



a) /(^) 



(2cj-l) 4(2cu + l) 



{2uj + 5) 3 MttGo 



+ c' 



(2tu + 3) 3(2cc; + 3)^ 
(2cj - 1) 4(2cj + 1) 



;iOcj + 13) 2 V 3c4 



a) ^(^)|r^ 



(2cj + 3) 3(2cu + 3)^ 



( 



AnGn 



A ^(^) r^ 



^ ^ -2/(2w+3) 
^0^ 



(Mcj + 19) V 3c4 

Thus again if we substitute a large-a; value, say, u ~ 10^ into eq.(29) above, we get, as 
f{u;),gicu)r^Oil), 



f^(r) 



100{km/s) X — 
ro 



^x -(1/106)^2 



(30) 



+ 



2n /G, 



-^ ^Ar, - +c'^ 






3 MttGo . -^NVr ^4 



10 V 3c4 ° 



-(1/106) 



3 V c^ / Vro/ iU \ dC V \ro. 

where we introduced the normalization factor tq ~ 10 (kpc) ~ 10^^ {cm) which represent 
a typical distance to the outer region of the dark halo where rotation curves begin to get 
flattened. We are now in a position to theoretically determine the upper bound on the value 
of the (positive) cosmological constant in our dark matter model, namely the BD theory 
in the presence of the cosmological constant. That is, in order for this result in eq.(30) to 
successfully describe the flattened rotation curve in the outer region of dark halo, r > tq, it 
should be imposed that 

OoA 



(A) 
(B) 



■^rl « (100A;m/s)^ 



OoA 



« 10" 
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as [GqK/c^] = l/s^ and [GqA/c^] = l/cm^. Then using Go = 6.67 x 10"^ (cmVc/s^) 
and c = 3 X 10^° (cm/s), the conditions {A) and {B) give A << 1.35 x 10~^ [erg/cm^). 
It thus is interesting to remark that this theoretical upper bound on the value of the 
cosmological constant is consistent with the observed value Aots — W^^ {erg/cm^) (or 



^obs/c^ — 10~^^ {g/cm?)) 171, [iSj. Of course an alternative interpretation is acceptable 
as well. That is, with the observed small value Kobs — 10~^ {erg /cm?), our model for dark 
matter turns out to be able to reproduce the flattened rotation curve in a successful manner. 
Therefore, it appears that our model, based on the Brans-Dicke theory in the presence of the 
cosmological constant, can successfully reproduce both the dark matter halo conflguration 
and the flattened rotation curves inside of it. This may indicate that our model could be 
one of the promising candidates for dark matter but other known evidences for dark matter 
need to be tested as well in this context of BD theory (possibly with cosmological constant) 
in order for it to be truly successful model of dark matter (and dark energy). 
Next, the equation of state in eq.(26) of this BD scalar fleld becomes, in the outer region of 
the galactic dark matter halo (i.e., at large-r). 



w 



2{uj - 1) cos2 e + sin^ o] + [(2a; - l){2uj + 3) + 4/3] (27rGo/c^) Ar^ sin^ 6 

(31) 



2{uj + 1) + [{2u - l){2uj + 3) - 4(2cj + l)/3] {2'kGq/c^) Ar^ sin^ 6 
[(2cj-l)(2cu + 3)+4/3] 



[{2iu - l){2iu + 3) - 4(2cj + l)/3] 

which is obviously negative due to the negative pressure (and still positive energy density) in 
this outer region. Moreover, for the large-a; value, i.e., uj ~ 10® for which the rotation curve 
gets flattened out, this equation of state at large-r further approaches w ~ — 1. Once again, 
this observation is particularly interesting as it appears to indicate that the BD scalar we 
are considering possesses dark energy-like negative pressure on larger, cosmological scales 
consistently with our previous study [41] that on the cosmological scale, the BD scalar field 
with A does exhibit the nature of dark energy possessing the negative pressure. 
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V. RIGOROUS RELATIVISTIC TREATMENT OF THE ROTATION VELOCITY 
AND THE BIPOLAR JETS 

A. The Rotation Velocity - a Fully Relativistic Derivation 

In the previous section, the rotation velocity of a test body orbiting in the background of 
the Schwarzschild-de Sitter-type spacetime in BD theory has been taken as the expression in 
the Newtonian hmit approximation. In the present section, the construction of the rotation 
velocity shall be promoted to a fully relativistic version in terms of the rigorous derivation 
of timelike geodesies in the Schwarzschild-de Sitter-type spacetime. 

The Schwarzschild-de Sitter-type spacetime in BD theory given in eq.(14) in the text is static 
and axisymmetric. Thus it possesses time-translational isometry generated by the timelike 
Killing vector ^^ = (d/dt)^ = 5f and rotational isometry generated by the axial Killing 
vector ip^ = {d/d(f))^ = 5'^. And the associated conserved quantities are energy E = E/m 
and angular momentum L = L/m (per unit rest mass) of a test body along the geodesic 

-E = M^e^ = ^„tw" = -(Asin2^) V--. T ) (rf^)' ^ ^ 

I = u,r = 9a^u^ = (a sin^ 9) ''^^'^^'^ r' sin^ 6 (^] (33) 

where u'^ = (dx^/dr) denotes the 4-velocity which is tangent to the geodesic. In addition, 
we have g^vU^u^ = —i^, namely 

-2/(2w+3) 



— K 



(a sin^ ^^ 
+ (Asm'' 9 



\.^A^^^Ac^(§)\r^.n.H(§)' 



-f-¥^f (ir-^m 



(34) 



where k = c^ for timelike geodesies (i.e., for massive bodies) and k, = for null geodesies 
(i.e., for light rays). Now, one may use eqs.(32) and (33) to eliminate (dt/dr) and {d(f)/dT) 
in terms of E and L and the result may be substituted into eq.(34) in order to obtain the 
"first integral" of the radial geodesic equation. Further, since gases or stars orbiting around 
the host galaxy center usually lie on the galactic (equatorial) plane, the case of equatorial 
geodesies at ^ = 7r/2 is of particular interest and the result is 

(^-j + Kjf{E,L;r) = 0, (35) 
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VeffiE,L;r) 



1 



2M SvrA 



+ ft;A-2/(2-+3) 



E' 



Therefore, now the problem of obtaining the timehke and null geodesies on the equatorial 
plane of the Schwarzschild-de Sitter-type spacetime in BD theory reduces to solving a prob- 
lem of non-relativistic, one-dimensional motion in an effective potential Veff{E,L]r). Here 
we consider the case of our particular interest, the "stable" circular orbit motions of massive 
objects characterized by the simultaneous conditions of 
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eff 



0, 



dK 



eff 



dr 







(36) 



as the ones first studied by Bardeen, Press and Teukolsky [19] for the case of Kerr spacetime. 
Then the conditions for the stable circular orbit motions in eq.(36) yield the required values 
of the specific energy and the specific angular momentum as 
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1- 



2M SttA 



^2^2 



3M/r) 
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'M 
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_ _^ ^2^-2/(2(^+3) 



SttA 
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167rA 



(37) 

A-2/(2-+3). (38) 



3 ; (2cj + 3) 

Note that it is the rotation angular velocity (per unit test mass) v^ rather than the specific 
angular momentum that we are after. Thus consider, again on the equatorial plane, the 
proper rotation velocity given by v"^ = {da^/dt) with ds^ = gudt^ + c/af (i.e., da^ denotes 
spatial line element) and dr = dO = 0, namely 



^-2/{2a;+3)^2_ 



{d(p/dT) 
(dt/dr) 



D/r^ 



(39) 



(1 - 2M/r - 87rAr2/3)-2E2/c4 
where we used eqs.(32) and (33). Further, substituting eqs.(37) and (38) into eq.(39) yields 



vl = c2A-2/(2-+3) X 



(40) 



(1 - 2M/r - 87rArV3) [{M/r - Mkr'^ /?,) - (2Z^(1 - M/r - 167rAr73)} 



(1 - 3M/r) + [{Mir - 87rAr2/3) - ji^^i^l - M/r - 167rAr2/3)} 

This is the fully relativistic expression for the rotation angular velocity (per unit test mass). 
Lastly, in the outer region of BD scalar halo, where r ~ 10 [kpc] ~ 10^^ (cm) >> GqM/c^ ~ 
10^*^ {cm) and hence A = [r'^{l - 87rArV3) - 2Mr]/rl ~ [r'^{l - 87rArV3)]/r2, this fully 
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relativistic expression for the rotation angular velocity reduces to 

(l-87rArV3) 



(1 - SvrArVS) - 
SttA 2 



(2w+3) 



;i - 167rAr73) 



(41) 



1- 



IGvrA 



r2 / 
- 1 



SttA 



-| -2/(2(^+3) 



3 (2cj + 3) V" 3 yj [ro 

Further, upon substituting a large-u; value, say, oj ~ 10^, one ends up with 

2 / o_ A \ ^ -2/(2a;+3) 



vl ~ C^ 



2 (2cu - 1) SttA 

(2cj + 3) ^ (2cj + 3) 3 



.^l^ 



SttA 



100(A;m/s) x <^ ^ 1 



'0 



SttG, 



+ 



,2 ^ 8^^,2 ) , _ 



3c4 



Kr' 



-(1/106)- 



^ttGo , ^ o 



(42) 



-Ar' 



-{1/10«) 



'0 V 3c4 

where as before the normalization factor ro ~ 10 {kpc) ^ 10^^ {cm) which represent a typical 
distance to the outer region of the dark halo has been introduced. As has been discussed 
in the text, therefore, in order for this result to successfully describe the flattened rotation 
curve in the outer region of dark halo, r > ro, it should be imposed that 

"'""""-d) « lo- 



se ■", 

which again yields A << 1.45 x 10~^ {erg/cm^). Note that this is essentially the same result 
as what we have gotten using the Newtonian mechanics version of the rotation velocity in 
section IV. To summarize, even if we employ the fully relativistic expression for the rotation 
velocity, the conclusions we have reached in the previous section essentially remain the same. 



B. Singularity along the Symmetry Axis: the Relativistic Bipolar Jet Interpreta- 
tion 

Earlier in section II, we mentioned in advance that a rigorous study of the true nature 
of the singularity along the symmetry axis reveals the fact that the bizzare singularity at 
^ = 0, vr of the Schwarzschild-de Sitter-type solution in BD gravity theory can account 
for the relativistic bipolar outflows (twin jets) extending from the central region of "active 
galactic nuclei (AGNs)". Thus in the present section, we shall demonstrate in an explicit 
manner that this is indeed the case. 
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In order to see if indeed the Schwarzschild-de Sitter-type spacetime solution in BD gravity 
theory allows for the occurrence of relativistic bipolar outflows, we particularly explore the 
geodesic motion of a test particle in the immediate vicinity of the symmetry axis, ^ = 0, vr. 
Namely, consider the geodesic motion at ^ ~ 5 << 1, (and hence sin^ ~ sin 5 = e << 1). 
Then the conserved quantities, i.e., the speciflc energy E = E/m and the speciflc angular 
momentum L = L/m of a test particle along the geodesic are now given by 



2 ,^-2/{2a;+3) / 2M SttA A ^ ( dt 



E^(.'A) -^.--^-■^r').'[^^^y»l, (43) 

Namely, the speciflc energy of a test particle is extremely large and its speciflc angular 
momentum is very small there for generic cj-value. And it is rather obvious to see that 

/ 2 \ -2/(2(^+3) 

mainly this can be attributed to the metric function factor ( A sin 6 ) (where sin 6 



is replaced by e ^ 1) which is responsible for the singular nature of the symmetry axis. 
Note flrst that the extremely small speciflc angular momentum near the symmetry axis is 
rather expected (like in the Einstein gravity context) as it is basically due to the negligible 
arm length from the symmetry axis. (Angular velocity (rf^/rfr) there, however, may not be 
so small.) The extremely large speciflc energy near the symmetry axis, however, is indeed 
something unexpected and thus is surprising as it does not happen in the Einstein gravity 
context (i.e., the a; ^ oo limit) where the symmetry axis is perfectly regular. Now, this 
study of the true nature of the singularity along the symmetry axis leads us to suspect that 
perhaps the bizzare singularity at ^ = 0, tt of the Schwarzschild-de Sitter-type solution in 
BD gravity theory can account for the relativistic bipolar outflows (twin jets) extending from 
the central region of active galactic nuclei (AGNs). That is, the curious singularity along 
the symmetry axis seems harmless, after all. Rather, it turns out to be a pleasant surprise 
as it can explain a long-known puzzle in observed features of galaxies. We now demonstrate 
that indeed this extremely large speciflc energy near the symmetry axis can be translated 
into the ultrarelativistic speed at which the test particle moves along the symmetry axis. 
To this end, we start with g^^u'^u'^ = —k, particularly near the poles, namely. 



/ 2a\"2/(2w+3) 



+ , 2^A2/(2-+3) 
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^ttA 






r 3 y ^ ^ 

where we used 5 « 1 and e << 1. We are now ready to study the first integral of the 
timehke (for massive particles, /t = c^) radial geodesic equation given by 
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Thus along the symmetry axis ^ = 0, tt, the j»roj)er ejection velocity would be given by 



(45) 



I dr^ 



(dr/dr) 
^" {dt/drf 



(46) 



2M SttA 



(2a;+l)/{2a;+3) 



[,2(^/^^)2]2/(2c.+3) 

/ 2M SttA 



where we used eq.(45) and A = [r^(l — 87rAr^/3) — 2Mr]/rQ in the second line and uj ~ 10® ~ 
large in the last line. Once again, it is interesting to note that if the value of the BD uj- 
parameter were not large enough (like o^ ~ 10® which has been assigned to reproduce the 
fiattened rotation curve in the galactic halos), the ejection speed of test particles along the 
symmetry axis (dr/dt) could encounter the danger of exceeding the speed of light due to 
the factor l/(e^)^/(^'^+'^) on the right hand side for (dr/dt). Lastly, for a typical galaxy 
with total mass M ~ 10^^ Mq and with the observed value of the cosmological constant 

1/9 

A ^ 10^^ (erg/cm^), we have 2GqM/c^ ~ 0.01 pc and (87rG'oA/3c'') ~ 4 Gpc and hence 
we end up with 



Vr 



C 1 



2M SttA 



1/2 



1 



(0.01 pc) 



1/2 



(47) 



r ^4 X 10^ pc) 

Observationally, it is well-known that the typical size of the active galatic nuclei is less than 
1 (pc) and the extent of typical galactic jets ranges from several (kpc) to a few (Mpc). Thus 
the ejection speed of test particles along the symmetry axis (i) just outside the AGN, i.e., for 
r ~ few (pc) is Vr — c[l — 0.01]^/^ ~ 0.995c and (ii) well-above the galactic plane, say, for r ~ 
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10 (kpc) is Vr — c[l — 10~^ — lO"^*^]^/^ < c. This completes the rigorous demonstration that 
the test particles move along the symmetry axis at nearly the speed of light. In addition, 
since the present analysis of the geodesic motion is valid only in the immediate vicinity of 
the symmetry axis, it indicates the collimation of the outflow as well consistently with the 
observation (see Fig. 2). 




FIG. 2: Centaurus A : An elliptical galaxy with an AGN and a jet. 



VI. CONCLUDING REMARKS 

In the present work, the pure (i.e., with no ordinary matter) BD gravity with or without 
the cosmological constant A has been demonstarted as a successful candidate for the theory 
of dark matter. To summarize, the mysterious flattened rotation curves observed for so 
long in the outer region of galactic halos, among others, have been the primary evidence 
for the existence of dark matter. Therefore, it has been demonstrated in this work that our 
model theory can successfully predict the emergence of dark matter halo-like configuration 
in terms of a self-gravitating static and nearly spherically-symmetric spacetime solution to 
the BD field equations and reproduce the flattened rotation curve in the outer region of 
this dark halo-like object in terms of the non-trivial energy density of the BD scalar field 
which was absent in the context of general relativity where the Newton's constant is strictly 
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a "constant" having no dynamics. As stated earlier, there has been a consensus that unless 
one modifies either the general relativity, the standard theory of gravity, or the standard 
model for particle physics, or both, one can never attain the satisfying understanding of 
the phenomena associated with dark matter and dark energy. Therefore, our dark matter 
model presented in this work can be viewed as an attempt to modify the gravity side alone 
in terms of the Brans-Dicke theory to achieve the goal. After all, the interesting lesson we 
learned from this study is the fact that the simplest extension of general relativity which 
involves relaxing Newton's constant to a dynamical field (i.e., the BD scalar field ^(x)) 
essentially makes all the differences. Thus to conclude, from this success of "BD scalar field 
as the dark matter" to account for the asymptotic fiattening of galaxy rotation curves while 
forming galactic dark matter halos plus the original spirit of the BD theory in which the 
BD scalar field is prescribed not to have direct interaction with ordinary matter fields (in 
order not to interfere with the great success of equivalence principle), we suggest that the 
Brans-Dicke theory of gravity is a very promising theory of dark matter. And this implies, 
if we emphasize it once again, that dark matter (and dark energy as well, see [^) might 
not be some kind of an unknown exotic "matter" , but instead the effect resulting from the 
space-time varying nature of the Newton's constant represented by the BD scalar field. Even 
further, this successful account of the phenomena associated with dark matter of the present 
universe via the BD gravity theory might be an indication that the truly relevant theory 
of classical gravity at the present epoch is not general relativity but its simplest extension, 
the Brans-Dicke theory with its generic parameter value u ~ 10^ fixed by the dark matter 
observation. This idea, however, should be taken with some caution as there are other 
phenomena observed thus far (than the fiattened rotation curves) which are suspected to 
be related to the effects of dark matter. These include the galaxy/cluster lensing of distant 
quasars. Besides, our model theory for dark matter presented in this work is not without 
a flaw. Earlier, from the rigorous study of the true nature of the singularity along the 
symmetry axis of this Schwarzschild-de Sitter-type spacetime, we suspected that perhaps 
this bizzare singularity at 6' = 0, vr can account for the relativistic bipolar outflows (twin 
jets) extending from the central region of active galactic nuclei (AGNs). The relativistic 
bipolar outflows, however, have been observed only for some particular types of galaxies 
such as AGNs and micro-quasars. That is, the jets do not seem to be a general feature 
of all types of galaxies. Therefore the Schwarzschild-de Sitter-type solution in BD gravity 
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theory employed in the present work appears to come with only a limited descriptive power 
for the observed phenomena associated with galaxies. At present in the absence of Birkhoff- 
type theorem guaranteeing the uniqueness of the Schwarzschild-type solution in BD gravity 
theory, however, it seems worth looking for an alternative solution and proceeding with it 
toward the issue addressed in the present work. 

Therefore more extensive and careful study needs to be done to test the Brans-Dicke theory 
as a truly successful model theory of dark matter and dark energy and we hope to report 
more along this line in the near future. 
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Appendix : Halos of BD scalar field : the case of Einstein frame 

Note that in the present work, we ignore the presence of other types of ordinary matters. 
Furthermore, we particularly consider the case when the cosmological constant term (which 
is a dimensionful parameter in the action) is neglected as well. Namely, we are dealing 
with the pure Brans-Dicke theory of gravity and we particularly restrict our interest to the 
original spirit of Brans and Dicke according to which the BD scalar field $ is prescribed 
to remain strictly massless without any self-interaction potential. Then in such a context, 
the study of possible role of BD scalar field as a dark matter component in Einstein frame, 
which is related to the Jordan frame (which we have been working in) by the Weyl rescaling 
(i.e., conformal transformation plus field redefinition), would be relevant as well. Thus we 
shall turn to this issue in this appendix. 
We now begin with the Weyl-rescaling given by 

g,, = n'ixyg,,, $ = My/"- (48) 

with n^x) = -^, ^o = (2t^ + 3) 
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where M^^ = l/Go denotes the Planck mass. Under this Weyl-rescahng, the action in the 
Jordan frame given earher in eq.(l) but without the A-term transforms to 



S 



J4 / - 

a xy g 



IGvr 



R - -r^'d.-^d,^ 



(49) 



Then extremizing this action with respect to the metric (^^j, and the redefined BD scalar 
field \E' yields the classical field equations given respectively by 
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V^^V,^ - -g,, [r'^'^^^Vp^ 



(50) 



The action and the classical field equations in this Einstein frame appear to take the forms 
of those of massless scalar field theory coupled minimally to Einstein gravity. A caution, 
however, needs to be exercised. That is, despite how it seems, this theory in the Einstein 
frame is really the (pure) BD gravity in disguise. Indeed, close inspection of eq.(48) above 
reveals that the Weyl rescaling becomes trivial only for cj — » oo, which is the Einstein gravity 
limit (with no scalar field), but for finite u values this theory is just the conformal rescaling 
of the (pure) BD gravity. As such, the solution to these field equations in the Einstein frame 
in eq.(50) should not be looked for independently without referring to that in the Jordan 
frame. Indeed, it would simply be the same Weyl rescaling of the solution constructed in 
the Jordan frame given in eq.(14), namely. 
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(51) 



We now attempt to address the same issues as we did earlier while working in the Jordan 
frame, i.e., whether the BD scalar field can cluster into dark matter halo-like objects and 
reproduce flattened rotation curves. Firstly, once again it is obvious that the BD scalar can 
cluster into halo-like configuration as it can be represented by this Brans-Dicke-Schwarzschild 
solution in the Einstein frame given in eq.(51). Secondly, we turn to the behavior of the 
rotation curve in the outer region of this halo-like configuration. To this end, again we 
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restore both the Newton's constant Gq and the speed of hght c to work in the CGS unit as 
before by multiplying the factor (c^/Gq) to the BD scalar field energy-momentum tensor in 
eq.(50) and by replacing M -^ GqM/c^ = M in the BDS solution in Einstein frame given 
in eq.(51). Next, by substituting the BDS solution in eq.(51) into the BD scalar energy- 
momentumm tensor in eq.(50) and then setting it equal to (15), again we can read off the 
energy density and the pressure components of the BD scalar field imperfect fluid as 



p = 


C^ 1 /. . 2 ^ -4/(2^+3) 


'{r-Mf cos^^ 
r2A ' r^sin^^ 


Pr = 


2.Go(2c. + 3)(^""') 


\r-Mf cos^^ 
r2 A r2 sin^ 9 


Pe = 


-Pr, P^ = -^P, 


Te = 


--'-^,j^M-'r'-'^^^^^- 



(52) 



Now, in order to study the behavior of the rotation curves in the outer region of 
dark halos of typical galaxies, we approximate the expression for the energy density 
of the BD scalar field at large-r, say, r >> GqM/c^ and it is p ~ c^(27rG'o(2c<J + 
3)) ""^(^^sin^ ^)^"'^ ( A sin^^) . Then the mass function is given by M{r) = 



Jq"^ d(f) J^^^ d9 Jq drJgrrgeeg(p(ppir,9) = {c'^/Go{2uj + 3)) f{5)r. Note here that we already 
have M{r) ~ r regardless of the value of BD c<j-parameter in the present case of Einstein 
frame and hence 

where this time f{5) = J^-'d9sm~^9 = 2 Jq^^ dx[l - x^]-^ ~ 0(1) with e, 5 « 1. This is 
already a flattened rotation curve, and if we further take, say, u ~ 10^ once again, we get 

v{r) ^ 100(A;m/s). (54) 

Therefore, regardless of whether we work in the Jordan or in the Einstein frame, the BD 

scalar fleld imperfect fluid always reproduces the flattened rotation curve. And from this 

observation, the BD scalar fleld appears to play the role of dark matter component in the 

galactic halos. 

Consider next, the equation of state of this BD scalar fleld imperfect fluid in this Einstein 

frame 

_ P _ [(r - M)Vr2 A] - [cos^ ^/r2 sin^ ^] 
^ ~ cV ~ [(r - M)2/r2A] + [cos2 9/r^ sin^ 9] ^ ^ 
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where P = Pr- Just like in the earher study in the Jordan frame, P = w{r,6)c'^p with 
w{r,6) ~ 0(1) meaning that this BD scalar fluid is still a barotropic fluid but with position- 
dependent coefficient w{r,6). And once again, although the BD scalar fleld is a candidate 
for dark matter, this equation of state indicates that in general, particularly in the vicinity 
of the individual galaxy, it is not quite a dust. In the outer halo region, i.e., at large-r, 
however, w ~ (as P ~ in the numerator) meaning that indeed the BD scalar fluid in 
the Einstein frame exhibits the property of dust matter everywhere including on the galactic 
(equatorial) plane regardless of the value of BD parameter uj. 

We now compare the results obtained earlier in the Jordan frame with those here in the 
Einstein frame. Firstly, in both frames, the BD scalar fleld imperfect fluid appears to 
play the role of dark matter component in a successful manner as it always reproduces the 
flattened rotation curve in the outer region of galactic halos. Secondly, the equations of 
state of the BD scalar fluid in both frames indicate that it is generally not quite a dust and 
this point implies that the BD scalar fleld appears to be an "exotic" type of dark matter 
possessing non-standard equation of state generally in the vicinity of individual galaxy. In 
the outer halo region, i.e., at large-r, however, it is interesting to note that lu ~ 0, namely a 
dust matter for large-cj parameter value and particularly on the galactic (equatorial) plane, 
where most of the gases or stars orbit around the galaxy center, in the case of Jordan frame 
whereas it is w ~ everywhere, not just on the galactic (equatorial) plane, and regardless 
of the value of uj in the case of Einstein frame. 

One might wonder how this somewhat delicate discrepancy in the nature of equation of 
state for the BD scalar fluid in the two conformal frames arises. We now attempt to provide 
one possible cause that may lead to this slight discrepancy. Indeed, close inspection reveals 
that it can be attributed to the non-standard form of the BD scalar fleld energy-momentum 
tensor in the Jordan frame. That is, in the expression for the BD scalar fleld energy- 
momentum tensor given in eq.(2), we can realize that it involves both the terms quadratic 
in flrst derivatives of the BD scalar fleld $ coming from its kinetic energy term in the action 
and those linear in second derivatives of $ coming from its non-minimal coupling term to 
curvature, ~ ^Jg^R. (Note that its standard expression for the ordinary scalar fleld (like 
that in the Einstein frame given in eq.(50) for the Weyl rescaled BD scalar fleld \1') involves 
only the terms quadratic in flrst derivatives of the scalar fleld.) And these anomalous terms 
linear in second (covariant) derivatives of $ seem to be the ones that make the difference as 
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they turn out to render the pressure of the BD scalar imperfect fluid negative at large-r and 
particularly for large-a; while zero on galactic equatorial plane in the Jordan frame. In both 
frames, however, the BD scalar field imperfect fluid appears to reproduce generic features 
of the dark matter in the outer region of dark halo of the individual galaxy. 
Note that in the other case where the cosmological constant is present, the issue of selecting 
from Jordan or Einstein frame that we just discussed becomes irrelevant to address. And 
it is because when a dimensionful parameter such as the cosmological constant (with mass 
dimension 4 in the Planck unit) is present in the theory action, the Weyl-rescaling of the 
action/field equations to go from the original Jordan frame, say, to Einstein frame does not 
leave the physics invariant. Therefore in section IV in the text, we just worked in the original 
Jordan frame. 
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